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^ i Abstract 

Annihilation rates of vector 1 charmonium and bottomonium 3 Si states V — > e + e~ 
' and V — ¥ 37, V — > jgg and V —¥ ig are estimated in the relativistic Salpeter method. In 

t-H , our calculations, special attention is paid to the relativistic correction, which is important 

and can not be ignored for excited 2S, 3S and higher excited states. We obtain T(J/ip — > 
^ 3 7 ) = 6.8 x 1(T 4 keV, T(ip(2S) -> 3 7 ) = 2.5 x 1CT 4 keV, T(ip(3S) ->■ 3 7 ) = 1.7 x 1CT 4 keV, 

Qh' r(T(15) 3 7 ) = 1.5 x 1CT 5 keV, r(T(2S') ->■ 3 7 ) = 5.7 x 1CT 6 keV, r(T(35) -> 3 7 ) = 

3.5 x 10- 6 keV and r(T(4S > ) -> 3 7 ) = 2.6 x lO" 6 keV. 

1 Introduction 

(N : 

> ; 

Annihilation decay of 1 S-wave heavy quarkonia has been extensively studied [1, 2, 3, 4, 5, 6, 7]. 
The interest in this study comes from several sources. First, the annihilation amplitudes are 

m ; 

related to the behavior of wave function, enabling an understanding of the formalism of inter- 
quark interactions. Further more, it can be a sensitive test of potential models [8]. Finally, the 
ratio of the decay widths, e.g. T(V — > e + e~)/T(V — > 3g), is sensitive to the running coupling 

> ■ 

constant a s (/x), where V is a heavy quarkonium vector state and /x is the scale (/i = m c for J/ip 
and [i = nib for T), and may provide very useful information for a s at the heavy quark mass 
scale [9, 10]. 

In our previous Letters [11], two-photon and two-gluon annihilation rates of J PC = ++ 
and 2 ++ cc and bb states are computed with the relativistic Salpeter method. Good agreement of 
the predictions with other theoretical calculations and the available experimental data is found. 
In the calculations, we found the relativistic corrections are large and not negligible, especially 
for high excited states, such as, the 25 and 3S states, because there are node structures in 
wave functions of 2S and 35 states, these cause large relativistic corrections even for heavy 
quarkonium like bottomonium. So in the theoretical studies concerning the highly excited 
states, a relativistic model is required. 
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The annihilation decays of the vector 1 states are different from the C even states. 
Basically there are two types of annihilation decay modes, which are V — > 7* — > l + l~ and 
V — > 37, jgg, 3g. These decay widths have been studied in non-relativistic limit and found to 
be proportional to the square of the wave function at the origin I ^(0)| 2 [1, 12]. However, the 
decay rates of many processes are subject to substantial relativistic corrections [10, 13, 14]. In 
this Letter, we will continue to study the annihilation decays of J PC = 1 cc and bb states 
with the relativistic Salpeter method. 

There are two sources of relativistic corrections [4, 11], one is the correction in relativistic 
kinematics which appears in the decay amplitudes through a well-defined form of relativistic 
wave function (i.e., not merely through the wave function at origin); the other relativistic cor- 
rection comes via the relativistic inter-quark dynamics, which requires a relativistic formalism 
to describe the interactions among quarks and relativistic formalism to consider the transition 
amplitude. To consider the relativistic corrections, we choose the Salpeter method [15], which 
is an instantaneous version of Bethe-Salpeter method [16]. For the equal-mass quarkonium, 
the non-instantaneous correction is very small [17]. For the annihilation amplitude, we choose 
Mandelstam formalism [18], which is well suited for the computation of relativistic transition 
amplitude with Bethe-Salpeter wave functions as input. 

In Section 2, we give theoretical details for the annihilation amplitude in Mandelstam for- 
malism and the corresponding wave function with a well-defined relativistic form. The decay 
width of V — > 7* — > e + e~ and V — > 37, jgg, 3g are formulated in this section. We will show 
the numerical results and give discussions in Section 3. 

2 Theoretical Details 

2.1 The V ->• e+e decay 

According to Mandelstam formalism [18], the transition amplitude of a quarkonium decaying 
into a electron and a positron (see figure 1) can be written as 

/d A q g^ v -> 

j<^^x{q)l^j^u ri {ki) lu v r2 (k 2 ), (1) 

where e q = 3 for charm quark and e q = — | for bottom quark; k\ and k 2 are the momenta 
of electron and positron respectively; M is the meson mass; %((/) is the Bethe-Salpeter wave 
function with the total momentum P and relative momentum q, related by 
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pi = aiP + q, ai = 



p 2 = a 2 P - q, a 2 = 



mi + ^2 
m 2 



mi + m 2 

where mi = m 2 is the constitute quark mass of charm or bottom. 

After performing the integration over q°, one reduce the expression, with the notation of 
Salpeter wave function *&(q) = i J ^xio)-, to 

/dq g^ v 
^^tr[^(^^j^u ri (ki)-f u v r2 (k 2 ). (2) 

We note that the form of the wave function is also important in the calculation, since the 
corrections of the relativistic kinetics come mainly through it. By analyzing the parity and 
charge conjugation, the general form of relativistic wave function of 1 _ state ( 1 for equal 
mass systems) can be written as [19] 



+ A?fe(q±) + {4±A ~ Q± ■ A)fr(q±) + ~ Pq± ■ A)fs(q±), (3) 

where P and are the momentum and polarization vector of the vector meson; q± = (0,q). 
The 8 wave functions are not independent due to the equations tp*r(q±) = p^(q±) = 0. For 
quarkonium states we get the constraints on the components of the wave functions [19]: 

x gj/3(g±) + M 2 / 5 (g ± ) 
h(q±) = ^ , MQ±) = 0, 

Mf 6 (q ± ) 

fs{Q±) = , fm±) = o. 

mi 

With these constraints, only four independent components f^, f^, f$ and fa are left. Namely 
^-M = <a • + |)/3^) + 1, ■ A^h(q,) 

+ {MA + 9± ■ ei-)h(q^) + [A P + ^ ' ^ ~ M^L],^). (4) 
mi mi mi 

These wave functions and the bound state mass M can be obtained by solving the full 

Salpeter equation with the constituent quark mass as input. We will not show the details of 

how to solve the full Salpeter equation, only give the final results. Interested readers can find 

the detail technique in Refs. [19, 20]. 

Defining the decay constant /y by 

f v Mel= (0|gi 7 ^ 2 |V,e) = y/l J 0-tr&(q)^}, (5) 



and with the Eq. (4) we can easily obtain 



Mo) Mo) 



(6) 



r)3 [ JO ™ 3M 2 

Summing over the polarizations of the final states and averaging over that of the initial state, 
neglecting the electron mass, it is easy to get the decay width 

Air 

r e+e - = -o?e\f Y lM. (7) 
2.2 V — y 37, V — y ^gg and V — y 3g decays 

With the notation and definition used in the previous subsection, the relativistic transition 
amplitude of a quarkonium decaying into three photons (see figure 2) can be written as 

T 37 = V3(iee q ) 3 f ^S^tTixiqM-n \ h~, \ h 

J (27r) 4 fo-p2-m pi-fti-m 

+all other permutations of 1,2,3]}, (8) 

where fci, &2, &3 an d ei, 62, €3 are the momenta and polarization vectors of three photons respec- 
tively. 

Since pio +P20 = M, we assume pio = P20 = M/2 as did in Ref. [11]. Having this assumption, 
we can perform the integration over qo to reduce the expression to 

d 3 q 



T 3l = V3(ee q fj A^tr{*(g)[/ 3 



#3— i>2 - rn f\- }t\ — m 
+all other permutations of 1,2,3]}, (9) 



where pi = {™,q),p 2 = {%,-<[)■ 
The decay width is given by 



M M 

The width of V — y jgg and V — y 3g are related to the three photon decay width by 

2 a 2 s 

r 79S = 3^| r3 7' ( n ) 

r 3<y = Tl^TT^S-y (12) 

54 cre° 

For gluonic decay V — >■ 3g, the trace of color generators gives tr [T a Ti,T c ] = \{d a h c + ifabc), so 
the expression of decay width contains two parts, one of which is proportional to the square of 
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symmetric constants and the other is proportional to the square of antisymmetric constants. The 
existence of antisymmetric term breaks the relation Eq.(12). However, our calculation shows 
that the antisymmetric term is sufficiently small compared to the symmetric term, so we can 
ignore it safely In non-relativistic limit the antisymmetric term vanishes exactly. 

3 Numerical Results and Discussions 

To solve the full Salpeter equation, we choose a phenomenological Cornell potential. There are 
some parameters in this potential including the constituent quark mass and one loop running 
coupling constant. The following best-fit values of input parameters were obtained by fitting 
the mass spectra for heavy quarkonium 1 states [21]: 

m c = 1.62 GeV, m b = 4.96 GeV (13) 

For cc system, we set Aqcd = 0.27 GeV. With this parameter set, we solve the full Salpeter 
equation and obtain the mass spectra shown in Table 1. To give numerical results , we need to 
fix the value of the renormalization scale fi in a s (fi). In the case of charmonium, we choose the 
charm quark mass m c as the energy scale and obtain the coupling constant a s {m c ) = 0.38 [11]. 

For bb system we set Aqcd = 0.20 GeV. With this parameter set, the coupling constant at 
the scale of bottom quark mass is a s (ra(,) = 0.23 [11]. The mass spectra are also shown in Table 
1. 

With the obtained wave function, Eq. (6) and Eq. (7), we calculate the decay width of 
V — > e + e~ for cc system. The results, with other theoretical predictions and experimental data 
from Particle Data Group, are shown in Table 2. Our results are larger than experimental 
data and consistent with the Beyer's [23] model version b results and Li's results [29]. The 
discrepancy between ours and experiment's may be due to the QCD corrections. We only have 
the leading order QCD correction 1 — [3] in hand, while the large factor ^ implies that 
high order QCD corrections can be still large and quite essential [30, 31], so we only show the 
results without QCD corrections. 

Decay widths of T(nS) — > e + e~ are shown in Table 3. All the results, with or without QCD 
corrections, are consistent with each other, with only small discrepancies. Since the small value 
of a s at the energy scale of bottom quark, corresponds to much smaller QCD corrections in 
bottomonium states than those in charmonium states, less discrepancies exit among the results 
of T(nS) decays than of Y>(nS) decays. 

The ratios of the high excited-state widths to the ground-state width r(nS)/r(lS) are free 
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from the QCD corrections and sensitive to wave functions. We show the ratios of leptonic decay 
widths in Table 4. Our theoretical values are comparable to the PDG data, except those in 
ip(3S) and ^(4S I ) states. It is can be seen from the table that the ratios, so do the decay widths, 
fall very slowly with successive radial excitations, which indicates that the relativistic corrections 
are large for high excited states. 

Decay widths T3 7 , T~ fgg and I^g of charmonia and bottomonia are calculated with Eqs. (9~12). 
The results and other theoretical estimates as well as experimental data are shown in Table 5 
and Table 6. The decay widths quoted from Ref. [9] and Ref. [30] are estimated based on exper- 
imental data. As in the e + e~ decays case, we only have the leading order QCD correction, e.g. 
1 — 12.6^ [3], in hand, and the large factor 12.6 implies that if we consider the QCD corrections, 
we need include high order QCD corrections not only the leading one. Besides, current leading 
order QCD factor is too sensitive to the value of a s , which make other contributions, such as 
relativistic corrections, unclear, so we only show the results without QCD corrections. For the 
same reason in the leptonic decays case, we show the ratios r3 9 (nS)/T3 9 (lS) in Table 7. 

A typical non-relativistic calculation gives B(J/ip — > 37) ~ 3 x 1CT 5 [28, 30], while our 
relativistic result #(J/Y> 37) ~ 0.73 x 10~ 5 (The total width of J/i/> is 93.2 ± 2.1 keV [22]) is 
much smaller than the non-relativistic one, but within the experimental error bar. This indicates 
that the relativistic corrections for charmonia 37 decays, so do the jgg, 3g decays, are large. 
This conclusion is also obtained by other authors, see Refs. [10, 13, 14]. Besides, our calculations 
show that the relativistic corrections for higher excited states are even lager than those for the 
ground state and lower excited states. 

One can see from Tables 5~7, that the decay widths, which are sensitive to the wave functions 
of corresponding states, fall very slowly from IS to 5S. We obtained the similar results as the 
cases of e + e" decays. This behavior is different from the non-relativistic models, where the values 
fall quickly from 15 to 5S. It shows that the relativistic corrections are large and important, 
especially for the higher excited states. It is believed that the relativistic corrections are small 
for bottomonium, however, we point that, this is true for ground state, but not exactly true for 
the excited states, especially for high excited states. 

In calculating the decay widths T^, we assume pio = P20 = M/2 . To show the effect of 
relaxing this assumption, we take p\o = 0.9 x M/2, P20 = 1.1 x M/2 and estimate the relative 
deviations of decay widths (T — To)/ro, which are shown in Table 8. We interchange the values 
of pio and P20, say, take pio = 1.1 x M/2, P20 = 0.9 x M/2, and find that the results are exactly 
the same as the unchanged case as expected. 



6 



In summary, by solving the relativistic full Salpeter equation with a well defined-form of wave 
function, we estimate annihilation decay rates of heavy quarkonium 1 ( 3 5i) states including 
V —> e + e~, V —> 37, V —> ^ygg and V — > 3g. We conclude that the relativistic correction 
and QCD correction in these annihilation decays play important roles, and high order QCD 
corrections are expected. 

Table 1: Mass spectra of the cc and bb 1 ( 3 <Si) states ( 3 Di states are not presented) in unit of 
MeV. The experimental data are taken from PDG [22]. 





Th(cc) 


Ex(cc) 


T(nS) 


Th(66) 


Ex(bb) 


J/V 


3096.9 


3096.916 


T(1S) 


9460.3 


9460.30 


V>(2S) 


3688.4 


3686.09 


T(2S) 


10024 


10023.26 


^(3S) 


4056.0 


4039 


T(3S) 


10371 


10355.2 


^(4S) 


4327.7 


4421 


T(4S) 


10635 


10579.4 


V>(5S) 


4543.3 




T(5S) 


10853 


10865 



Table 2: Decay width r(^(nS) —> e + e ) in unit of keV. The results marked by f do not cover 
the contributions of QCD corrections. 



i/}(nS) 


J/1> 


^(2S) 


V>(3S) 


V>(4S) 


V>(5S) 


Ours^ 


10.9 


5.22 


3.49 


2.61 


2.07 


Beyer Ver.a^ [23] 


5.33 


2.31 


1.59 


1.14 




Beyer Ver.b T [23] 


11.2 


4.06 


2.74 


2.06 




EQt [24] 


8.00 


3.67 








VPBK^ [25] 


5.469 


2.140 


0.796 


0.288 




GVGV [26] 


2.94 


1.22 


0.76 


0.43 


0.27 


IS [27] 


6.72 ± 0.49 


2.66 ±0.19 








SYEF [28] 


3.93 


1.78 


1.11 






LOT [29] 


11.8 


4.29 


2.53 


1.73 


1.25 


PDG [22] 


5.55 ±0.14 ±0.02 


2.38 ± 0.04 


0.86 ± 0.07 


0.58 ±0.07 






Figure 1: Leptonic decay diagram of quarkonium. 
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Table 3: Decay width T(T(nS) — > e + e ) in unit of keV. The results marked by f do not cover 
the contributions of QCD corrections. 



TfnS) 












Ours"!" 


1.47 


0.736 


0.530 


0.425 


0.359 


Beyer Ver.a^ [23] 


1.24 


0.51 


0.35 


0.28 




Beyer Ver.b 1 " [23] 


1.41 


0.56 


0.36 


0.30 




EQT [24] 


1.71 


0.76 


0.55 






VPBKt[25] 


1.320 


0.628 


0.263 


0.104 


0.0404 


GVGV [26] 


0.98 


0.41 


0.27 


0.20 


0.16 


IS [27] 


1.45 ± 0.07 


0.52 ±0.02 


0.35 ±0.02 






Gonzalez [32] 


1.7 


0.61 


0.39 


0.27 


0.21 


PDG [22] 


1.340 ±0.018 


0.612 ±0.011 


0.443 ± 0.008 


0.272 ±0.029 


0.31 ±0.07 



Table 4: The ratios of the high excited-state widths to the ground-state width r(nS)/T(lS) for 
decay T(^(nS) -> e+e") and T(T(nS) -> e+e"). 



r^(nS))/r(^(is)) 


r(2S)/r(is) 


r(3S)/r(is) 


r(4S)/r(is) 


r(5S)/r(is) 


Ours 


0.48 


0.32 


0.24 


0.19 


PDG [22] 


0.43 


0.15 


0.10 




r(T(nS))/r(T(lS)) 


r(2S)/r(is) 


r(3S)/r(is) 


r(4S)/r(is) 


r(5S)/r(is) 


Ours 


0.50 


0.36 


0.29 


0.24 


PDG [22] 


0.46 


0.33 


0.20 


0.23 



Table 5: Decay widths of ip(nS) — > 37, ^ygg, 3g in unit of keV. The data marked by * is quoted 
from Ref. [35]. The results marked by f do not cover the contributions of QCD corrections. 



Decay 


Ours^ 


GP [8] 


ML [33] 


POP 1 " [34] 


SYEF [28] 


Voloshin [30] 


r(J/V -> 3g) 


101 


176 


80 ±40 


63.72 




61.5 ±3.1 


r(^(2S) -> 3g) 


36.6 


78.4 




20.49 




45.3 ±9.3 


r(V(3S) -> 3g) 


24.7 






11.92 






r(^(4S) -> 3<?) 


19.8 






8.08 






r(^(5S) ->• 3 5 ) 


16.7 












r(j/V> -> 755) 


6.18 




7.5 ±3 






7.46 ±2.80 


r(^(2S) -> 755) 


2.25 










3.04 


r(V(3S) -> 755) 


1.52 












r(^(4S) -> jgg) 


1.22 












r(^(5S) -> 755 ) 


1.03 












r(j/V> -»• 3 7 ) 


0.68 x 10~ 3 








0.56 x 10~ 3 


*(1.12±0.47) x 10~ 3 


r(^(2S) -> 37) 


0.25 x 10~ 3 












r(^(3S) 37) 


0.17 x 10~ 3 












r(^(4S) -> 37) 


0.13 x 10~ 3 












r(^(5S) 37) 


0.11 x 10~ 3 
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Table 6: Decay widths of T(nS) — > 37, ^ygg, 3g in unit of keV. The results marked by f do not 
cover the contributions of QCD corrections. 



Decay 


Ours^ 


nrt roi 
KjL 1 [oj 


Mij [00J 


KlVixtK, [yj 


TVYYl _v 1n\ 

1 1 (lz>) — > og) 


oz.o 


AA 1 


OO _|_ C 


A 9 Q _|_ 1 9 


r(T(2Sl -> 3ol 


12.0 


22.5 








7 47 


1 6 9 










12.1 








4.41 








1 ( I (loj -> jgg) 


U.OZD 




n 4- n o 
u.y ± U.Z 


1 9n 


T(T(2S) -> 755) 


0.304 








T(T(3S) -> 755) 


0.190 








T(T(4S) -> 755) 


0.140 








r(T(5S)^ 7 ^) 


0.112 








r(T(is) -»• 37) 


0.15 x 10~ 4 








T(T(2S) -»• 3 7 ) 


0.57 x 10~ 5 








r(T(3S) -> 3 7 ) 


0.35 x 10~ 5 








T(T(4S) -> 3 7 ) 


0.26 x 10~ 5 








T(T(5S) -»• 3 7 ) 


0.21 x 10~ 5 









Table 7: The ratios of the high excited-state widths to the ground-state width r(nS)/T(lS) for 
decay r(^(nS) -> 3g) and T(T(nS) -> 3g). 



r(v>(nS))/r(v(is)) 


r(2S)/r(is) 


r(3S)/r(is) 


r(4S)/r(is) 


r(5S)/r(is) 


Ours 


0.36 


0.24 


0.20 


0.17 


GI [8] 


0.45 








PCP [34] 


0.32 


0.19 


0.13 




Voloshin [30] 


0.74 








r(T(nS))/r(T(lS)) 


r(2S)/r(is) 


r(3S)/r(is) 


r(4S)/r(is) 


r(5S)/r(is) 


Ours 


0.37 


0.23 


0.17 


0.14 


GI [8] 


0.51 


0.38 


0.27 





Table 8: Relative deviations of decay width with assumptions p±o = 0.9 x M/2, P20 = 
1.1 x M/2 from that with p w = p 2 o = M/2. 



nS 


IS 


2S 


3S 


4S 


5S 


cc 


-2.0% 


-2.9% 


-4.7% 


-7.0% 


-8.2% 


bb 


-1.8% 


-2.0% 


-2.2% 


-2.4% 


-2.6% 
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Figure 2: 37, jgg, 3g annihilation diagrams of quarkonium. We do not show other diagrams 
with all the possible permutations of photons and gluons. 
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